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ABSTRACT 

We have identified an obvious and strong large scale structure at redshift 
z 5.75 in a wide (31'x33') field, narrowband survey of the Chandra Deep Field 
South region. This structure is traced by 17 candidate Lya emitters, among 
which 12 are found in an 823nm filter (corresponding to Lya at 2; = 5.77 ± 0.03) 
and 5 in an 815nm image = 5.70 ± 0.03). The Lya emitters in both redshift 
bins are concentrated in one quadrant of the field. The Hubble Ultra Deep Field, 
Chandra Deep Field South, and GOODS-South fields all lie near the edge of this 
overdensity region. Our results are consistent with reports of an overdensity in 
the UDF region at 2; ~ 5.9. This structure is the highest redshift overdensity 
found so far. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: 
high-redshift — large-scale structure of universe 



1. Introduction 

The study of the large-scale clustering of high redshift galaxies is essential to understand 
the formation and evolution of galaxies. To date, many possible detection of large-scale clus- 
tering of high redshift galaxies have been reported (e.g., Steidel et al. 1998, 2000; Venemans 
et al. 2002, 2004; Miley et al. 2004; Palunas et al. 2004; Ouchi et al. 2001, 2003, 2004; 
Foucaud et al. 2003; and references therein). 
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The LycK emission line searches have been very effective at finding high redshift galaxies, 
especially at 2; > 5 where most of the spectroscopically confirmed galaxies show Lja emission. 
These include our Large Area Lyman Alpha (LALA) survey (e.g., Rhoads et al. 2000) and 
other recent searches (Cowie & Hu 1998; Hu et al 1998, 2002, 2004; Kudritzki et al 2000; 
Fynbo, MoUer, & Thomson 2001; Pentericci et al. 2000; Stiavelh et al. 2001; Ouchi et al. 
2003; Fujita et al. 2003; Shimasaku et al. 2003; Kodaira et al. 2003, Ajiki et al. 2004; Kurk 
et al. 2004; Tran et al. 2004, and references therein). Large-scale clustering of high redshift 
Lyo; emitters has also been reported in some of these surveys (Venemans et al. 2002, 2004; 
Ouchi et al. 2001, 2003, 2004). 

In this letter we report the detection of a large-scale structure of Lya emitters at 2; ~ 5.7 
in Chandra Deep Field South (CDF-S). The candidate Lya emitters at 2; ~ 5.7 were selected 
using deep narrowband images, following the selection criteria developed in the Large Area 
Lyman Alpha survey (Rhoads et al. 2000, 2001, 2003, 2004). Throughout this paper, we 
assume a cosmology with Hq — 71 km s~^ Mpc"-*^, JIm = 0.27 and = 0.73 (cf. Spergel et 
al. 2003). 

2. Deep Narrow Band Imaging 

We imaged the area around Chandra Deep Field South using two narrowband filters 
with central wavelength Ac of 815 and 823 nm. The full width at half maximum (FWHM) 
transmissions of the two filters is 7.5 nm, with peak throughputs of ~ 90%. These were new 
filters of high quality. In particular, their throughput does not vary substantially across the 
field of view (unlike the similar but older filters used by Rhoads & Malhotra 2001). The 
deep narrowband images were obtained using the Mosaic 11 CCD imager at the CTIO 4m 
V. M. Blanco telescope, on 2003 Dec. 22-30 (UT). Chip 3 of the 8 SITe 2Kx4K CCDs 
was turned off during the observations due to high amplifier noise. Imaging data reduction 
followed the methods described in Rhoads et al. (2000, 2004). To summarize, in addition 
to the standard CCD data reduction steps (overscan subtraction, bias frame subtraction, 
and flat- fielding) , we also removed electronic crosstalk between Mosaic chip pairs sharing 
readout electronics, and the residual, large-scale imperfections in the sky flatness using a 
smoothed supersky fiat derived from the science data. Cosmic rays were rejected in each 
exposure using the algorithm of Rhoads (2000). Astrometry of USN0-A2 catalog stars was 
used to adjust the world coordinate systems of individual frames. Satellite trails were fiagged 
manually for exclusion from the final stacks. Weights for image stacking were determined 
using the "ATTWEIGHT" algorithm (Fischer & Kochanski 1994), and the task "mscstack" 
in the MSCRED package (Valdes 1998) was used to stack the individual exposures. 
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The total integrated exposure time of the final stacked images are 34.7 ks for the 815 
nm filter and 36.0 ks for the 823 nm filter. The final stacked images have seeing of 0.9", 
while the seeing of 823 nm band image is slightly better (0.02") than that of 815 nm band 
image. Wc also used deep B, V, I broad band images of ESO Imaging Survey (EIS) in the 
CDF-S, obtained using the Wide Field Imager (WFI) of the ESO-MPG 2.2 meter telescope 
at La Silla. The data were download from ESO archive, and registered to the narrow band 
images for candidate selection. The overlap between the narrowband stacks and the EIS 
images was 31 x 33 arcminutes in size with a gap of 10.5 x 13 arcminutes due to the missing 
chip. The total sohd angle searched for Lya galaxies was thus 887 square arcminutes. The 
corresponding volume over the full redshift range, 5.67 < z < 5.80, is 3.06 x lO^Mpc^ A 
summary of the images is given in Table 1. 



3. Candidate z— 5.7 Lya Emitters 

Source detections were performed using SExtractor version 2.2.2 (Bertin & Arnouts 
1996) on the narrowband images, and measured their colors using SExtractor's two-image 
mode to measure the EIS B,V,I band photometry for the narrowband detected sources. All 
fluxes are measured in 2.43" (9 pixel) diameter aperture. 

Candidate z ~ 5.7 Lya emitters were selected following the same criteria that have 
proven highly successful in the LALA fields (Rhoads & Malhotra 2001; Rhoads et al. 2003, 
2004). These are: 1) Narrowband detection > 5(t; 2) A narrowband excess of > 0.75 mag- 
nitude, so that > 50% of the narrowband fiux comes from an emission line; 3) Significance 
of the narrowband excess > 4(t; and 4) < 2a detection in filters blueward of the expected 
Lyman break location (5, V band in this letter). The success rate of the criteria has been 
spectroscopically confirmed to be > 70% (Rhoads et al. 2003, Dawson et al. 2004, Dawson 
et al. 2005 in prep). 

In order to make full use of our narrow band images, which are deeper than the / band 
image, wc ran our selection using the weighted mean fiux from I and 815 nm band as the 
underlying continuum for 823 nm selection, and the weighted mean fiux from / and 823 nm 
band as the continuum for 815 nm selection. A total of 17 candidates are selected, 5 in the 
815 nm band, and 12 in the 823 nm band. In Fig. 1 we plot the underlying continuum and 
narrow band fluxes for all detected sources and candidate z ~ 5.7 Lya emitters selected. 
Our selection criteria are also plotted in the flgure: the vertical line stands for 5 a detection 
in narrow band; the solid inclined line represents a narrow excess of 0.75 magnitude and the 
dashed line corresponds to the requirement of 4cr significance of the narrowband excess. 
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4. Lcirge scale structure 

Among the 17 candidate z ~ 5.7 Ly« emitters we detected in CDF-S, 12 of them are 
selected in the 823 nm band image, 5 of them in the 815 nm band image. In Fig. 2 we 
plot the sky distribution of the 17 candidate z ~ 5.7 Lya emitters in CDF-S. An obvious 
structure of the distribution can be seen in the figure: almost all the 17 sources are located 
in the lower left half of the field. 

Before studying the clustering properties of the 17 Lya emitters, we consider whether 
the structure could be an observational artifact. Such a spurious "structure" in the sky 
distribution of the selected Lya emitters could be due to two reasons: It might be due 
to a difference in narrow-band sensitivity across the field. It might also due to the sky 
distribution of very bright sources in the field, which could affect our source detection and 
candidate selection. 

To check these possibilities, we first divide the 823 nm band image into four quarters 
and plot the source counts for the four parts in Fig. 3. We can see that the source counts 
in four regions are consistent within the Poisson counting uncertainties. In region a where 
we identified the most Ly« candidates, the overall source density is among the lowest. This 
indicates that the detection efficiency in region a is not higher than in other 3 regions. We 
found a similar result in the 815 nm band image. 

Second, we added artificial point sources with narrowband magnitude of 24.3 (the mean 
value of the 12 Lya emitters in the 823 nm band image) at random locations in the narrow- 
band image. We ran exactly the same procedures of source detection and candidate selection 
on the derived narrowband images. The recovered artificial sources passing our criteria show 
no evidence of strong clustering. Inserting artificial sources with different magnitudes (the 
highest and the lowest magnitudes of the 12 sources) yields the same results. Likewise, no 
sign of clustering is found among artificial sources inserted and recovered from the 815 nm 
band image. A two-dimensional Kolmogorov- Smirnoff test shows that the sky distribution 
of the 12 Lya emitters at 823 nm is different from the corresponding simulation at the 99.4% 
confidence level. Including the 5 Lya emitters from the 815 nm band would slightly increase 
the significance to 99.6%. The sky distribution of the 5 Lya emitters in 815 nm band alone is 
different from the simulated sources at the level of 98%. We checked the limiting magnitudes 
of the 823 nm, /, B and V band images over the whole field, and found no inhomogeneity. 
The photometric zero points of the narrow and broad band images are also found to be 
accurate within 0.1 mag over the whole field. We point out that there is no evidence of 
significant spatial variations in bandpass of the narrow band filters either. All the above 
checks indicate the clustering of the Lya candidates is not artificial. 
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We see clear evidence of large scale clustering of the 12 Lya emitters in the 823 nm band. 
The same trend is also seen in the 815 nm band, but with lower confidence level because the 
number of Lya emitters in 815 nm band is only one third of that in 823 nm band. In Fig. 4 
we compare the source counts from the 815 nm and 823 nm band images. We can see that 
the detection efficiencies arc consistent with each other above the 5 a detection limit, which 
is 24.45 mag for the 815 nm band image, and ~ 0.2 mag deeper for the 823 nm band. Note 
there are 2 Lya emitters in 823 nm band with narrow band magnitudes fainter than 24.45. 
Removing the 2 sources, the space density of Lya emitters in 823 nm band is still twice of that 
in 815 nm band. The corresponding limiting Lya fluxes are 2.00 & 1.66 X 10-^'^ ergs /cm'^/s 
for 815 nm and 823 nm bands respectively. Based the luminosity function of z ~ 5.7 Lya 
emitters given by Malhotra & Rhoads (2004), we predict the numbers of Lya emitters in the 
two narrow band images to be 4.5 & 7.5. This corresponds to an ovcrdcnsity factor of 1.6 
for 823 nm band candidates. In the densest quadrant of the survey, the overdensity factor 
reaches 3-4. 



5. Conclusions 

The structure we have found in the Lya galaxy distribution at z ~ 5.8 is the most 
distant large scale feature yet reported in the spatial distribution of galaxies. The space 
density of these galaxies varies by a factor of two between our two narrowband filters, and 
the projected surface densities span a substantially larger range (up to 4x the expected 
mean) . 

The Hubble Ultra-Deep Field (HUDF; Beckwith et al 2005 in prep) lies near the edge of 
this large scale galaxy sheet, and the effects of this ovcrdcnsity arc also apparent there. An 
overdensity at 2; ~ 5.9 near the HUDF was first suggested by Stanway et al (2004) on the 
basis of 3 redshifts in that region. This has now been confirmed in the redshift distribution 
of 23 i-dropout Lyman break galaxies in this field (Malhotra et al 2005), as measured using 
slitless spectroscopy from the GRAPES project (Pirzkal et al 2004). The redshift extent of 
the feature within the UDF is larger than the coverage of our narrowband filters, suggesting 
that the structure is more complex than a single sheet. 

We would like to thank Bahram Mobasher for providing the cafibrated photometric zero 
points for the EIS images. The work of JW was supported by the CAS "Bai Ren" project 
in USTC and the CXC grant G03-4148X. 
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Table 1. Optical images 



Band 


Telescope 


Exposure Time (ks) 


mAB(lim)" 


B 


ESO 2.2m 


24.5 


27.13 


V 


ESO 2.2m 


9.2 


25.97 


I 


ESO 2.2m 


26.9 


24.62 


815nm 


CTIO 4m 


34.7 


25.03 


823nm 


CTIO 4m 


36.0 


25.14 



^ Limiting magnitudes (AB) for a 3 cr detection on a 
2.0" diameter aperture. 
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Fig. 1. — Underlying continuum fluxes vs narrow band fluxes for SExtractor detected sources 
in the two narrow band images. Our candidate z 5.7 Lya emitters (5 in 815 nm band, and 
12 in 823 nm band) are marked as big dots. Lines show our selection criteria: the vertical 
line stands for 5 a detection in narrow band; the solid inclined line represents a narrow excess 
of 0.75 magnitude {Fluxnarrow/ Flux continuum = 2) and the dashed line corresponds to the 
requirement of 4(T signiflcance of the narrowband excess which is FluXnarrow - FluXcontinuum 
= 4x a/ {Err'^^„.^,^^, + Err'^^^^^^^^^) . Note there is one candidate in the left panel which does 
not pass the three threshold lines in the figure. This is because that the lines are plotted 
using the mean uncertainties of continuum and narrow band fiuxcs, while these candidates 
are selected based on their individual flux uncertainties output by SExtractor, which can 
be 10-20% smaller/larger than the mean uncertainties. The negative underlying continuum 
fluxes in the flgure are just for display purpose. While selecting candidates, these negative 
fluxes were promoted to zero to make our selection conservative. 
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Fig. 2. — Sky distribution of the 17 photometrically selected z ~ 5.7 Ly« emitters in CDF-S. 
The shaded region in the plot is due to the Chip 3 of the Mosaic II camera which was turned 
off during the observations. North is up and east is to the left. Solid and open circles are 
candidates selected from 823 nm and 815 nm band respectively. See text for details. 
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Fig. 3. — The source number counts in the 823 nm band image, which was divided into 
4 quarters. The numbers of detected sources in each quarter agrees well with each other, 
showing that the clustering of our 2; = 5.7 candidates is not due to the different detection 
efficiency over the field. (The counts in quadrant c are adjusted for the reduction in solid 
angle due to the missing chip.) 
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Fig. 4. — Comparing the source number counts in the 815 nm and the 823 nm band images. 
The vertical hnes stand for the 5 a detection in the narrow bands as one of our selection 
criteria. 



